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Abstract 

We present an analysis of a partially directed walk model of a polymer which at one 
end is tethered to a sticky surface and at the other end is subjected to a pulling force at 
fixed angle away from the point of tethering. Using the kernel method, we derive the full 
generating function for this model in two and three dimensions and obtain the respective 



<4-H 
O phase diagrams. 

We observe adsorbed and desorbed phases with a thermodynamic phase transition in 

between. In the absence of a pulling force this model has a second-order thermal desorption 

transition which merely gets shifted by the presence of a lateral pulling force. On the other 

hand, if the pulling force contains a non-zero vertical component this transition becomes 

first-order. 

O Strikingly, we find that if the angle between the pulling force and the surface is beneath 

a critical value, a sufficiently strong force will induce polymer adsorption, no matter how 

►^ large the temperature of the system. 

Our findings are similar in two and three dimensions, an additional feature in three 

dimensions being the occurrence of a reentrance transition at constant pulling force for small 

temperature, which has been observed previously for this model in the presence of pure 

£ — . vertical pulling. Interestingly, the reentrance phenomenon vanishes under certain pulling 

angles, with details depending on how the three-dimensional polymer is modeled. 



• !-i 1 Introduction 

X 

There is ongoing interest in the study of linear polymers and their conformal entropy using 
directed and partially directed walk models, both from theoretical [U [2] and experimental [3] 
perspectives. These models are of particular interest because it is often possible to obtain exact 
solutions for their generating functions, partition functions and free energy. Thus a precise 
understanding of the phase transitions that the models can exhibit may be obtained, which can 
then be used to guide and interpret laboratory experiments on physical polymers [3]. 

In appropriate physical settings, for example when the polymer undergoes surface adsorp- 
tion, good qualitative agreement has been shown between the phase behavior of directed walk 
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models and that of self avoiding walk (SAW) models |3], the latter of which exhibit more of the 
configurational possibilities available to real polymers, but are less tractable as mathematical 
models. Under the presence of a pulling force, this is intuitive, because the application of such 
a force to a physical polymer tends to straighten it, making it behave more like a directed or 
partially directed path [5]. 

In this paper we consider a two-dimensional partially directed walk model for a polymer 
which at one end is tethered to a sticky surface and at the other end is subjected to a pulling 
force in a fixed direction away from the point of tethering. We then extend this model into 
three dimensions in two different ways. The first extension is 'minimally three-dimensional', in 
that its projection onto the surface is fully directed, and as a consequence walks adsorbed in 
the surface are fully directed. This extension has previously been considered in [6] for a purely 
vertical pulling force. The second three-dimensional extension we consider is different from the 
first in that walks adsorbed in the surface are partially directed. To our knowledge, this model 
has not been considered previously. 

Another novelty in our model is in the variable angle with which the pulling force may be 
applied to our tethered polymers. In similar previously published literature the pulling force was 
restricted to be vertical El El El El OH E] ; a limitation which does not apply to the technology 
of optical tweezers which is used by experimentalists in 'pulling-force' experiments carried out 
on real polymers. Changing the pulling direction has been considered previously in anisotropic 
self-interacting models [12l [131 E] • 

Our use of partially directed walks, rather than simpler directed models such as Dyck, 
Ballot and Motzkin paths, is for the purpose of allowing the extra degree of freedom required to 
enable polymers to extend horizontally when pulled upon with a force containing a horizontal 
component. 

2 Partially Directed Walks 

Let W = Zx Z>o be the square integer lattice in the upper half plane. Let p = Voeov±ei...eLVL 
be a path of length L consisting of an alternating sequence of vertices and edges such that each 
vertex Vi £U and each edge e^ := (uj_i, v{) is of one of the following three types: 

(0,1) a vertical step up 
Vi-i = ^ (1,0) a horizontal step across (2-1) 

(0,-1) a vertical step down 

We impose three additional constraints on walks in our model. The first is tethering to the 
origin, imposed by the condition vo = (0, 0). The second is that paths are self-avoiding, that is, 
an up step may not be immediately followed by a down step, or vice versa. We call such paths 
tethered partially directed walks. The third constraint is that paths must always start and end 
on a horizontal step. This latter constraint makes negligible difference to the physics, and is 
introduced for computational convenience in finding the generating function. 



Geometrically relevant parameters of a walk are: 



N 
M 
K 
R 
H 



the number of horizontal edges in the walk, 

the total number of vertical edges in the walk, 

the number of horizontal edges of the walk lying in the surface, 

the number of vertical edges of the walk leaving the surface, 

the final height of the walk. 



We define the weight of a partially directed walk to be 



x N y M K K uV. 



(2.2a) 
(2.2b) 
(2.2c) 
(2.2d) 
(2.2e) 



(2.3) 



An example of a partially directed walk and its weight is illustrated in Figure [TJ Also illustrated 
is a pulling- force at angle 6 away from the iV-axis. 



HA 



Partially directed walk ending horizontally 
tethered at origin and 
beginning horizontally 



Pulling force 




Figure 1: A partially directed walk with weight x n y < k w 3 // 3 . For clarity in the diagram, the 
weights on the edges leaving the surface are indicated visually by dots on their initial vertices. 



Note that the conjugate pair of variables u and R, which keeps track of vertical departures 
from the surface, is introduced for later convenience in our consideration of our second three 
dimensional model. For the purposes of the two-dimensional and first three-dimensional model, 
these variables are unnecessary, and setting u = 1 in the following calculations simplifies them 
considerably. 

Also note that setting uj = k corresponds to considering a model in which vertices of the walk 
that are in the surface are weighted. We have considered both edge and vertex weighted models, 
and have found negligible differences. For simplicity we only present results for the edge- weighted 
models in this paper. Results for vertex-weighted models follow mutatis mutandis. 



3 Exact Solution of the Generating Function 



We define a generating function 
G(x,y,K,u;n) 



N,M,K,R,H>0 



CN,M,K,R,HX N y M K K U R jJL H 



(3.1) 



in terms of the path weight function given in Equation ( 2.3 ) , where c/v m K R H counts the number 



of different configurations with given parameters N, M, K, R, and H. A functional relation 
for the generating function may be obtained by considering paths ending at a fixed height. 
Suppressing the first four variables, we abbreviate G{(j) := G(x,y,n,uj; (J,). A combinatorial 
decomposition of the set of all walks with respect to generation by addition of 'hooks' of arbitrary 
height and ending in a horizontal step leads to the following functional equation. 



G(ji) 



h.v 



horizontal step at height (3.2a) 




+ G(y)( K -l)x 
+ G(0)(w-1 



JJf 



l~W 



horizontal step at height > 
vertical steps up, then horizontal 

vertical steps down, then horizontal 

removes walks descending below the surface 

plus contact weights for edges in the surface 
plus contact weights for edges leaving the surface 



(3.2b) 
(3.2c) 

(3.2d) 

(3.2e) 

(3.2f) 
(3-2g) 



The derivation of the functional equation is summarized in Lines (3.2a)-(3.2g). Since our 



paths are required to start and end in a horizontal step, the shortest possible path consists of a 



single horizontal step, and is accounted for by Line (3.2a). Lines (3.2b)-(3.2d) account for longer 



paths created by concatenating paths ending at height H, as encoded by the factor of G(fi), with 



paths containing a single horizontal step and any number of vertical steps. Line (3.2d) entails 



some over-counting, as it includes paths which contain vertices below the upper half plane, U. 



This over-counting is compensated for by Line (3.2e), which removes from the count any paths 
containing edges strictly below the iV-axis. 

So far paths which touch but don't drop below the surface have been included without a 



weight k on those edges in the path that lie in the surface. Line (3.2f ) corrects for this. Similarly. 



Line (3.2g) corrects for the missing weight u on vertical edges whose initial vertex lies in the 



surface. 

The functional equation thus obtained is amenable to the so-called kernel method [15]. First 
we simplify the functional equation somewhat by observing that G(0) and G{y) are simply 
related by 

G(0) = Kx(l + G(y)), (3.3) 



due to the fact that any walk ending with a horizontal step in the surface must either be a 
single step walk or have been obtained by dropping H steps down from height H and adding a 
horizontal step in the surface. 



Using Equation 3.3, the functional equation is expressed in Equation (3.4) below in a form 



suitable for solution by means of the kernel method, 



kernel 



KX I 1 + (UJ — l)x 



yv> 



i-yv 



l 



i-y/yu 



K — (uj — 1)KX 



Ufi 



i -yn 



G(y). (3.4) 



The coefficient of G(fi) is called the kernel. Setting the kernel equal to zero gives a quadratic 
equation for /j,: 



yv 



(1 - x + y 2 + xy 2 )fi + y = 



(3.5) 



Of the two roots, one is physically meaningful - we call it \jl v and denote the non-physical 
root by [i np . The correct root may be identified by expanding the power series and choosing the 
one with correct asymptotic behavior for G as y — > 0. Equivalently, the correct choice may be 
identified as the one which is consistent with the combinatorial interpretation of G{y). Using 
this criterion we choose the root whose series expansion contains no negative exponents, and 
has only non-negative coefficients in that expansion. Choosing the correct root fi = fi p gives 



G(y) 



K I 1 + (OJ — l)x- 



Wp_ 

-yv P 



l 



K — (CO — 1)HX- 



Wp_ 

-yv P 



(3.6) 



Thence, back-substituting G{y) into Equation (3.4) gives the full general solution for G(fi); 
explicitly: 



G(x,y,K,u;n) 



1 + (uj- l)x 



W 



KX- 



i-yn 



1 + (w - l)x 



y^ P 
1 -yv P 



l 



i-y/> 



K — ( UJ — 1)KX- 



!Jf' 



1-W 



1 



i - y/np 



K — [UJ — 1)KX- 



Wp_ 

i -y^p 



where 



x(l - y 2 ) 



l-y{li+l)+y 2 



Mp 



(1 — x + y 2 + xy 2 ) — y / — 4y 2 + (1 — x + y 2 + xy 2 ) 2 

2y 



(3.7) 
(3.8) 



4 The Path-Length Generating Functions 



In this section we introduce the changes to the generating function (3.7) that are needed for 
our three different models. We aim to write the generating function such that we can start to 
analyze the finite-step partition functions Zl, for walks of length L, which are coefficients of t L 
in an expansion G = X^l>o t Z l , where t is a new variable that is conjugate to the path length 
L. 

4.1 The 2D model 

The finite-step partition function Zl, for walks of length L in two dimensions, is the coefficient 
of t in the following expansion of the path-length generating function 

G(x = At, y = t, k, co = 1; fi) = ^ t L Z L (X, fi, k) (4.1) 

L>0 

under the substitutions 

(x, y, k, w; fj) -<— (At, t, k, 1; fi) , (4-2) 

so that t is conjugate to the path length, A is conjugate to the horizontal position and fi is 
conjugate to the vertical position. In this way we only give a contact weight to edges in the 



surface. In Section 5.1 we analyze the singularities of Equation (4.1 ) to determine the asymptotic 
growth of Zl in two dimensions. 

An alternative model would be obtained by weighting all vertices in the surface, leading 
to the consideration of G(\t, t, k, k; fj). We have completed both analyses and have found no 



significant difference, hence we only present the analysis of the case defined by Equation (4.1). 



4.2 The first 3D model 

Our first three-dimensional model is obtained by a substitution which replaces any step in the 
^-direction with a step in either the N\ or iV2 direction. The projection of the resulting walk 
onto the N\ — N2 plane is fully directed. This model is equivalent to the bi-colored walk model 
considered in [Sj. This corresponds to a substitution in the generating function variables given 
by 

(x,y,K,uj;iJ,) -«- (Ait + A 2 t, t, k, 1; fj,) (4.3) 

so that t is conjugate to path length, Ai is conjugate to the iVi-coordinate, A2 is conjugate 
to the ^-coordinate and \i is conjugate to the //-coordinate. A three-dimensional partially 
directed walk of our first kind is illustrated in Figure [2] and the asymptotic growth of Zl in 



three dimensions is analyzed in Section 5.4 below 



4.3 The second 3D model 

Our second three-dimensional model is obtained by a substitution which has the following effect. 
Given a two-dimensional partially directed walk in the N — H-plane, we insert arbitrarily long 
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Figure 2: A tethered, partially directed walk in 3 dimensions, on a sticky surface with contact 
weights k within the N\ — N 2 surface and a pulling-force applied at fixed angle <j> away from 
the A^-axis and angle 8 away from the N\N 2 surface. This walk ends at position (iVi, N 2 , H) = 
(8,4,0) and has weight A?A| t 21 K 5 /j°. 



sequences of steps perpendicular to the plane containing the original walk, in positions prior to 
each step in the original walk, while observing self- avoidance in the three-dimensional walk thus 
created. In this three-dimensional context we rename the original A-axis with new name, N\, 
and create a new axis, A2, perpendicular to both the N\ and H axes. The projection of the 
resulting walk onto the N\ — N 2 plane is no longer self-avoiding; likewise with the N 2 — H plane. 
The projection of the resulting walk onto the N\ — H plane is partially directed. The set of valid 
walks in our second 3D model is a strict superset of the set of valid walks in our first 3D model. 

In order to construct these three-dimensional walks we introduce a TL-coloring of the two- 
dimensional partially directed walks that we have already defined, by associating with each edge 
a color drawn from the set Z. These Z- colored walks are somewhat analogous to the two-colored 
walks in [9|, in which every horizontal edge was assigned one of two colors. In our context 
we assign all possible edge-colorings to all possible walks. Now, given a walk with a fixed 
coloring, we insert before an edge of color k, a sequence of |fc| edges in the positive or negative 
^-direction, as given by the sign of /c. In this way, every coloring of a two-dimensional walk 
produces a unique three-dimensional walk. An example of a Z-colored walk is given in Figure [3j 
The corresponding three-dimensional walk is illustrated in Figure [4} 

The associated substitution rules for the weights can be derived as follows. In the bulk (that 
is, for edges whose initial vertex does not lie in the surface), a fc-colored vertical edge of weight 
y gets replaced by 1 + \k\ edges of combined weight £ 1+ l fc l, and summing over all colors leads to 
the substitution 



y <- t (... + t 2 + t + 1 + t + t 2 + ...)= t 



1 + t 
1-t 



(4.4) 



and a fc-colored horizontal edge of weight x gets replaced by 1 + \k\ edges of combined weight 
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Figure 3: (Color online) A Z-colored walk with colors —2, — 1, 0, 1, 2 occuring, where —2 is shown 
in blue, —1 is shown in green, is (implicitly) black, 1 is red and 2 is orange. 
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Pulling force 




Figure 4: A tethered, partially directed walk in 3 dimensions, on a sticky surface with contact 
weights n within the A r i — N2 surface and a pulling-force applied at fixed angle 9 within the 
N\ — H plane. This walk ends at position (Ni,N2,H) = (7,4,0) and has weight \ 7 t sl K 5 /j,°. 



At 1+ I fc l, leading to 

x *- At-^- . (4.5) 

1-t v ; 

In the surface (that is, for edges both of whose vertices lie in the surface), a /c-colored edge of 

weight kx gets replaced by 1 + \k\ edges of combined weight \Kt(Kty k ', and summing over all 

colors leads to the substitution 

kx <— Ak£ . (4.6) 

I -Kb ' 

Finally, a A;-colored edge that leaves the surface, having weight ojy, gets replaced by 1 + 1 k\ edges 
of combined weight t(nty ', and summing over all colors leads to the substitution 

uy^t—. (4.7) 

The generating function for this three-dimensional model (or, equivalently, for the Z-colored 
walks), is therefore given by the following substitution 

/ 1 + t 1+t (l-t)(l + Kt) (l~t)(l + Kt) \ /A , 

(^>^) <- { xt T^t> 't^< K (iT7) (T^t)' (i+T) jr^ry n ■ (48) 

5 Singularity Analysis and Phase Diagrams 

Phase transitions in the two and three dimensional models, as path length tends to infinity, 
occur when singularities of the generating function closest to the origin coincide. 

We calculate the singularities of the generating function by considering the discriminant 
of the quadratic equation (G — G P )(G — G np ) = G 2 — (G p + G np )G + G p G np satisfied by the 
generating function, where G p and G np are respectively the physical and non-physical solution, 



corresponding to substituting roots n p and // np of Equation (3.5) into Equation (3.7). 
The discriminant, in original variables, is 



where 



and 



A = AiA 2 (5.1) 



x 2 (y - 1) (y + 1) (y + xy - x + 1) (y + xy - 1 + x) 

(-M + y + yp - y v -x^y 1 + x\i) 



k 2 uj 2 (Kxyn(y - n - yu> + ^uj) - ny(l - yfi + fi 2 ) - /x(l - k - yn)) 
(k 2 x 2 (1 — ujy 2 )(l — uj) + K0jy 2 (l + x — kx) — kx(2 — uj — kuj) — kuj + 1) 

5.1 The 2D model: singularity analysis 



In two dimensions, with path-length generating variables given by substitutions (4.2) into Equa- 



tion (5.1), the discriminant of the quadratic for G becomes 

K 2 X 2 t 2 (fit - l) 2 ( K t - (« - l)^) 2 (t -l)(t + 1) (At 2 - (A - l)t + 1) (At 2 + (A + l)t - 1) 



A 



1 \ \ 2 / \ 2 



(jl(k - 1) ( At 3 + t 2 - ( A + \i + - j t + 1 ) I K\t 3 1 2 - nXt + 1 



(5.4) 



We find that the relevant singularities of G in t are given by roots of 

At 2 + (A + l)t-l 



At 3 + t 2 



n\r 



(5.5a) 
(5.5b) 

„-i - < 5 ' 5c ' 

The first gives rise to algebraic square-root singularities, whereas the latter two give rise to 
simple poles. 

When [i = 1 in (5.5b), it factorizes as 



A + /U+- ) t + 1 
-t 2 -nXt + 1 



(i-l)(At 2 + (A + l)t-l) 



(5.6) 



and thus contains (5.5a) as a factor. Thus for fi = 1 the algebraic singularities coincide with 



poles given by the roots of (5.5b). This coincidence implies the occurrence of divergent square- 



root singularities (as opposed to convergent square-root singularities when fj, ^ 1). When // > 1, 



the pole given by the smallest root in t of (5.5b) dominates the algebraic singularities. 



The simple poles arising from the smallest roots in t of (5.5b) and (5.5c) coincide when 

,2\ 



A 



K/J,(k — 1 — K[i 



K 2 H 2 ] 



(5.7) 



(«- !)[(«- I) 2 

which defines a surface in the space defined by A, \x and k, and is shown in Figure [5] This surface 
separates the space into two regions and a boundary between them. In one region, the smallest 



root in t of (5.5b) is the singularity closest to the origin. In the other region, the smallest root 



in t of (5.5c) is the singularity closest to the origin. On the boundary surface, both of these 



singularities coincide. To interpret the meaning further, it is convenient to change to physical 
variables. 



5.2 The 2D model: physical variables 

The relevant physical variables for our system are the temperature T of the statistical mechanical 
ensemble of polymers, and the force F by which the polymer is pulled at a fixed angle with 
respect to the horizontal axis. This leads to horizontal and vertical force components F x = 
FcosO and F y = Fsin8, respectively. The energy E of a polymer under the influence of a 
pulling force is given by 

E = KJ - NF X - HF y (5.8) 

where KJ = (number of contacts) x (energy per contact), NF X = (horizontal distance) x 
(horizontal force), and HF y = (vertical distance) x (vertical force). 
The Boltzmann weight of such a configuration is then given by 



e -E/kT _ e -{KJ-NF x -HF y )/kT 



(5.9) 
= -1, 



which after scaling, such that the Boltzmann constant k = 1 and the attractive energy J 
becomes 

e -E/kT = K K X N ^ (5 1Q) 
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Figure 5: The surface in phase-space given by Equation (5.7), relating A, fi and k 



where 






e T 



f 



F sin I 

e t 



(5.11) 



The critical surface in terms of the physical variables, whose equation is determined by 
substituting Equations ( |5.11[ ) into Equation (5.7), is illustrated in Figure [6} 

The surface in Figure[6]is indicative of a phase transition between an adsorbed and a desorbed 



state for the polymer being modeled. This interpretation will be confirmed in Section 5.3 where 
we plot the fraction of the polymer which is in contact with the surface, for various pulling angles 
between 0° and 90°. As will be seen, that part of phase space which is behind the critical surface 
(as pictured in Figure pi) corresponds to a state of adsorption, in which a positive fraction of the 
polymer is in contact with the surface. That part of phase space which is in front of the critical 
surface as pictured corresponds to a state of desorption. 



5.2.1 General observations 

Some features of the desorption transition become clearer when considering slices through the 
surface at constant pulling angle, leading to temperature- force desorption curves as shown in 
Figure [7j 

When no force is applied, the polymer is adsorbed at low temperature and becomes desorbed 
as the temperature increases past a critical value. This thermal desorption occurs at 



T = l/log(l + V2/2) « 1.87 



(5.12) 
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0.0 



Figure 6: The 2D model: critical surface in physical variables T, F and 6. Slices through this 
surface at constant angle 9 are shown in Figure (TJ 



F 2 




Figure 7: The 2D model: temperature-force desorption transitions at different pulling angles 
from purely vertical to purely horizontal pulling in increments of 9° . The left-most curve corre- 
sponds to vertical pulling at 90°, and the right-most curve corresponds to horizontal pulling at 
0°. 
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as can be computed by setting force equal to zero, i.e. letting \i = A = 1 in Equation (5.7). 
Naturally, in the context of zero force, there is no dependence on the angle. When a force is 
applied and pulling is in a purely vertical direction, the force favors desorption, as expected. 
This may be seen by tracing upwards either through the front of the phase surface as pictured 
in Figure [6j or equivalently through the left-most curve within the temperature-force plots in 
Figure [7J (The area to the left of a curve in the temperature- force plot corresponds to a state of 
adsorption; and the area to the right to desorption.) As is also intuitive, pulling with a purely 
horizontal angle favors adsorption, since horizontally stretched polymers will favor adsorption. 
This may be seen by considering either the back of the phase surface as pictured, or the right- 
most curve in the temperature- force plot. 

5.2.2 Critical angles 

Perhaps surprisingly, we see that there is a non-zero critical angle below which pulling will never 
induce desorption no matter how great the force, and will eventually induce adsorption. To 
determine this value, we compute the angle for which the temperature- force curve has vertical 
slope at F = 0. The force-temperature plots in Figure [7] show this value to be about 27°; and 
in fact the critical angle is 

9 = tan" 1 (1/2) « 26.6° . (5.13) 

At zero temperature, the transition to desorption occurs when 

' (5-14) 



sin v — cos a 
Hence, once < 45°, no zero-temperature desorption occurs. 

5.3 The order parameter 

To confirm that there is indeed a transition between desorbed and adsorbed phases, we consider 
the fraction 

C= lim -(K) (5.15) 

L— >OQ Lr 

of the polymer that is in contact with the surface, which is an order parameter for the transition. 

In the desorbed state, we expect C = 0, whereas it will be positive in the adsorbed state. For 

a polymer lying entirely in the surface C = 1. A second-order phase transition between an 

adsorbed and a desorbed state will be reflected in a smooth change of C, whereas a first-order 

phase transition will be indicated by a jump of C from zero to a non-zero value. 

We calculate C from 

dlogt c 
C = -rr , (5.16) 

OlogK 

i.e. from the change of the critical fugacity t c with respect to the interaction weight k . In 
Figure [8] we show the temperature- force dependence of C at various fixed pulling angles ranging 
from vertical to horizontal in decrements of 18°. For zero force, we see that C changes continu- 
ously with temperature, indicative of a second-order phase transition. As shown in the bottom 
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Figure 8: The order parameter C, the fraction of polymer adsorbed onto the surface, as a function 
of temperature T and force F, for pulling angles 9 = 90°, 72°, 54°, 36°, 18° and 0° from top left 
to bottom right. 
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right diagram of Figure [8j the transition remains second-order when a horizontal pulling force 
is applied, and only the location of the transition shifts. On the other hand, as shown by all 
the other diagrams of Figure [8j the transition becomes first-order as soon as there is a non-zero 
vertical component of the pulling force. Note that once adsorbed, the value of C is independent 
of the pulling force F and angle 0. 

Thus, Figure [8] confirms that there is indeed a desorption phase transition, as indicated 
above. 



5.4 The first 3D model: physical variables 



According to the substitution (4.3), we now have variables Ai conjugate to horizontal variable 



N% and A2 conjugate to the other horizontal variable iVJz. The singularity analysis that was 



carried out for two dimensions in Section 5.1 is nearly unchanged, since it amounts to replacing 



A with Ai + A2 in Equation (5.4) 



The physical variables for three dimensions are now given by substitutions 



C'T 



Ai 



F cos cos <p 



M 



F cos sin d> 



/' 



(5.17) 



where, as before, T is temperature, F is force and 6 is pulling angle upwards away from the 
base surface. The extra variable is the component of pulling-angle measured away from the 
iVi-axis in the iViiVVplane. 



F 2 




F 2 




Figure 9: The first 3D model: shown are temperature-force desorption transitions for horizontal 
pulling angle = 0° (left) and 4> = 45° (right). In each diagram the curves correspond to 
different vertical pulling angles 8 from purely vertical to purely horizontal pulling, in increments 
of 9°. The left-most curve corresponds to vertical pulling at 90°, and the right-most curve 
corresponds to horizontal pulling at 0°. 



As for the two-dimensional case discussed above in Subsection 5.2, we consider temperature- 
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force desorption curves at constant vertical pulling angles 0. In Figure [9] we show the resulting 
diagrams for two different horizontal pulling angles = 0° and eft = 45°. The pictures are 
qualitatively similar to Figure [6] for two dimensions, with a few interesting differences. 

The phenomenon of reentrance is observable in both diagrams of Figure [9j as is seen by the 
shape of the temperature- force curves for steep angles. As may be observed in either diagram 
of the figure, when pulling vertically with a force F slightly larger than one, the polymer is 
desorbed at sufficiently high temperature, gets adsorbed upon decreasing the temperature, but 
then desorbs again when the temperature is decreased even further towards zero. 

Notably, there is a difference between the two diagrams, in that in the left diagram (</> = 0°) 
the reentrance phenomenon becomes weaker and disappears upon decreasing 9, whereas in the 
right diagram (<p = 45°) the reentrance phenomenon persists up until a critical value of 6, at 
which zero-temperature desorption disappears completely, is reached. 

Reentrance has been observed before in related models, see e.g. |llj . An explanation of the 
occurrence of reentrance is provided by a zero-temperature entropy argument, which shows that 
the critical force changes for small temperature T linearly as 

F^l + TS (5.18) 

where S is the the configurational entropy available to the walks adsorbed onto the two- 
dimensional surface. 

In this three-dimensional model walks adsorbed onto the surface are fully directed. When 
pulling at a horizontal angle cj) = 45°, all 2 L possible adsorbed configurations of walks of length 
L are equally likely, and hence the configurational entropy log 2 is positive. This effect is clearly 
visible in the temperature-force diagram via the positive slope of the curves for low temperature. 

When pulling at a horizontal angle cp 7^ 45°, there is a preferred direction for the stretched 
walk, hence the configurational entropy is zero. This effect is visible in the temperature-force 
diagram for eft = 0, where a non-vertical pulling angle immediately leads to a horizontal slope 
of the curves for low temperature. We argue that this effect is ultimately responsible for the 
disappearance of the reentrance phenomenon. 

As for the two-dimensional model, we can calculate several quantities exactly. For example, 
the critical value of thermal desorption is given by 

T = 1/ log(7/8 + y/rf/8) « 3.03 . (5.19) 

For (j) = 45° the non-zero critical angle below which pulling cannot induce desorption and will 
eventually induce adsorption is given by 

9 = tan~ 1 ((l + a/17/17)V2/4) « 24° , (5.20) 

zero-temperature desorption occurs when 

F= 2 — , (5.21) 

2sinfl-\/2cos6> 

and the critical angle below which no zero-temperature desorption occurs is 

6» = tan~ 1 (\/2/2) « 35° . (5.22) 
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5.5 The second 3D model: singularity structure and physical variables 



After the substitution (4.8), we obtain variables as in the two-dimensional model, i.e. we only 
vary the vertical pulling angle 9 (analogous to letting <fi = in the first 3D model). The 



singularity analysis that was carried out for two dimensions in Section 5.1 needs to be done 



carefully, as the substitution (4.8) potentially changes the phase diagram. The resulting algebraic 



equations become rather large and cumbersome, for example the equation corresponding to 



the surface given by Equation (5.7) for the two-dimensional model now becomes an algebraic 



equation in A, k, and /i which involves 932 monomials. We are therefore restricted to performing 
a numerical analysis. 

The physical variables are 






A 



F cos £ 

e t 



/' 



F sin 

e t 



(5.23) 



where, as before, T is temperature, F is force and 9 is pulling angle upwards away from the base 
surface. 




Figure 10: The second 3D model: shown are temperature-force desorption transitions for differ- 
ent vertical pulling angles 9 from purely vertical to purely horizontal pulling, in increments of 9°. 
The left-most curve corresponds to vertical pulling at 90°, and the right-most curve corresponds 
to horizontal pulling at 0°. 

We again consider temperature- force desorption curves at constant pulling angles 9 away from 
the N\ — N2 plane. In Figure [TU| we show the resulting diagram. The picture is qualitatively very 
similar to the diagram for pulling at = 0° in Figure [9] for the first three-dimensional model. 

The phenomenon of reentrance is again observable in Figure [9| Similarly to what was 
observed in the temperature- force diagram of the first three-dimensional model for (p = 0, 
a non-vertical pulling angle 9 immediately leads to a horizontal slope of the curves for low 
temperature and ultimately to the disappearance of the reentrance phenomenon for sufficiently 
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shallow pulling. We argue that this is again due to the fact that, with pulling as described, a 
stretched polymer adsorbed onto the surface at T = has zero configurational entropy. 

6 Conclusion 

Naively, one might expect that pulling sufficiently hard on a polymer tethered to a horizontal 
surface, at any angle containing a non-zero vertical component away from that surface, would 
eventually affect desorption. We have shown that in our models this is not the case, and that 
pulling with any positive angle below a certain critical value will induce adsorption. With 
our study of three-dimensional models, we have also shown how the reentrance phenomenon, 
which is present in force-induced polymer desorption, is affected by changing the entropy of 
the adsorbed polymer due to pulling with a horizontal force component. Interestingly, in semi- 
flexible polymers, which are modeled on a lattice by weighting bends, hence making the polymer 
stiffer, reentrance is affected in a somewhat similar manner [9J. It would be interesting to extend 
our study by including the effect of stiffness. 
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